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Abstract 

Capsazepine, an antagonist of capsaicin, is discovered by the structure and activity relationship. In previous studies it has 
been found that capsazepine has potency for immunomodulation and anti-inflammatory activity and emerging as a 
favourable target in quest for efficacious and safe anti-inflammatory drug. Thus, a 2D quantitative structural activity 
relationship (QSAR) model against target tumor necrosis factor-a (TNF-a) was developed using multiple linear regression 
method (IVILR) with good internal prediction (r^ = 0.8779) and external prediction (r^pred = 0.5865) using Discovery Studio 
v3.5 (Accelrys, USA). The predicted activity was further validated by in vitro experiment. Capsazepine was tested in 
lipopolysaccharide (LPS) induced inflammation in peritoneal mouse macrophages. Anti-inflammatory profile of capsazepine 
was assessed by its potency to inhibit the production of inflammatory mediator TNF-a. The in vitro experiment indicated 
that capsazepine is an efficient anti-inflammatory agent. Since, the developed QSAR model showed significant correlations 
between chemical structure and anti-inflammatory activity, it was successfully applied in the screening of forty-four virtual 
derivatives of capsazepine, which finally afforded six potent derivatives, CPZ-29, CPZ-30, CPZ-33, CPZ-34, CPZ-35 and CPZ- 
36. To gain more insights into the molecular mechanism of action of capsazepine and its derivatives, molecular docking and 
in silico absorption, distribution, metabolism, excretion and toxicity (ADMET) studies were performed. The results of QSAR, 
molecular docking, in silico ADIVIET screening and in vitro experimental studies provide guideline and mechanistic scope for 
the identification of more potent anti-inflammatory & immunomodulatory drug. 
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Introduction 

Capsicum species commonly known as cliUlies, relislied as an 
important spice in the culinary art of the world and is known for its 
medicinal effect since the dawn of the human civilization. The 
medicinal property of 'hot pepper' has been attributed to the 
presence of capsaicin, a pungent principal ingredient produced as 
a secondary metabolite. Chemically known as 8-methyl-N-vanillyl- 
6-nonenamide. Capsaicin and its related compounds, collectively 
referred as 'capsaicinoids' or VaniUoids', which bind specifically to 
transient receptor cation channel subfamily V (TRPV), that carry 
sensation of pain and responds naturally to noxious stimuli like 
high temperature and acidic pH [1]. Prolonged exposure causes 
nociceptor terminals to become insensitive to capsaicin, as well to 
other noxious stimuli [2]. Hyper stimulation of TRPVl by 
capsaicin has an analgesic effect, since it leads to long-term 
desensitization of the sensory neurons. The clinical uses of TRPVl 
agonist like capsaicin, are limited due to side effects of a burning 
sensation, irritation and neurotoxicity [3]. On the other hand, 
blocking of the pain-signalling pathway with a TRPV 1 antagonist 



capsazepine represents a promising strategy for the development 
of novel analgesics with potentially fewer side effects [4] . Several 
non-neuronal effects of capsaicin have also been reported viz., 
induction of apoptosis in transformed cells [5], stimulation of 
prostaglandin formation leading to inhibition of gastric lesion [6], 
antibacterial activity [7], inhibition of cardiac excitability [8] and 
platelet aggregation [9]. Capsazepine is a known analog of 
capsaicin, discovered as a result of structure-activity relationship 
(SAR) studies [10]. Capsazepine induced similar action as 
capsaicin and resiiiiferatoxin (RTX) and exhibits even twofold 
more potent inhibition of expression of iNOS gene in LPS- 
stimulated murine macrophages through inactivation of NF-kB 
[11,12]. NF-kB is a protein complex that control transcription of 
DNA and it is involved in cellular responses to stimuli such as 
stress, cytokines, free radicals, ultraviolet irradiations, oxidized 
low-density lipoprotein, and microbial antigens. NF-kB regulation 
of immune response and inflammation, cell lineage development, 
cell apoptosis, cell cycle progression and oncogenesis in response to 
stimuli have been shown to regulate the expression of several genes 
(bcl-2, bcl-xl), cellular inhibitor of apoptosis protein, tumor 



PLOS ONE I www.plosone.org 



1 



July 2014 I Volunne 9 | Issue 7 | e100797 



QSAR Studies of Capsazepine Derivatives for Immunomodulatory Activity 



80 



70 



60 




0 10 20 30 40 50 60 



Experimental IC50 \iM 

Figure 1. Graphical plot between experimental and predicted activities (IC50 ^M) of the training and test set compounds. (A) Training 
data set (blue dots), (B) Test data set (orange dots). 
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necrosis factor signalling pathway-related regulatory factor, 
cyclooxygenase-2 (COX-2), matrix metalloprotein peptide-9 
(MMP-9) and inducible nitric oxide synthase (iNOS) and those 
for cell cycle regulatory components involved in tumorigenesis 
[13-18]. Therefore, in this work we have investigated the chemo 
preventive potential of capsazepine and its derivatives against pro- 
inflammatory mediator TNF-ot through QSAR, in vitro activity 
evaluation and molecular docking studies, to understand the 
mechanism of action of vanilloids against inflammation and 
immunomodulation related to cancer. QSAR modelling also 
furnished the activity dependent structural descriptors and predicts 
the effective dose of other derivatives, thereby suggesting the 
possible toxicity range. Drugabihty of hit compounds was 
evaluated by using Lipinski's 'Rule of Five' and ADMET analysis 
through bioavailability filters. 

Materials and Method 

Dataset preparation 

A total of 146 known TNF-Qt inhibitors were collected from 
various pubhshed literatures based on its structural diversity and 
activity coverage. The activity data for all compounds were taken 
from difiFerent scientific groups in terms of inhibitory activity (IC50 
|J,M) [19-24]. 124 compounds out of 146, were selected as a 
training set based on following criteria to produce a good 
quantitative QSAR model: by covering a wide activity range of 
compounds and by including the most active, moderate and less 
active inhibitors (Table SI in File SI). The biological activity for 
TNF-Q( inhibitors were ranging between 0.09 to 50 |J,M. The 
remaining 22 compounds were used as a test set to validate the 
generated model (Table S2 in File SI). 



Energy minimization 

The structural drawing and geometry cleaning of the training 
set compounds were performed through, ChemBioOffice suite 
Ultra vl2.0 (2010) software (CambridgeSoft Corp., UK). The 
compounds then subjected to energy minimization by using 
Discovery Studio v3.5 software (Accelrys Inc., USA) by applying 
CHARMm forcefield applicable to most of the small molecules. It 
adds several properties to the compounds including: initial 
potential energy, RMS gradient, CHARMm energy and mmimi- 
zation criteria. 

Chemical descriptors calculation 

Molecular descriptors were calculated for each compounds 
using "Calculate Molecular Properties" module of the Discovery 
Studio v3.5 (Accelrys Inc., USA). These descriptors include 2D 
parameters [e.g., AlogP, molecular weight, number of aromatic 
ring, number of H-acceptors, number of H-donors, number of 
rings, number of rotatable bonds, molecular fraction polar surface 
area) and 3D (Dipole and Jurs descriptor). 

Quantitative structure activity relationship (QSAR) model 
development 

The set of energy optimised 146 compounds with calculated 
molecular properties were used for QSAR model development 
using create QSAR model module in Discovery Studio v3.5. 
Firstly all compounds were prepared for QSAR, and then the 
biological activities were specified as dependent property. Com- 
pounds were randomly divided into training (124 compounds) and 
test (22 compounds) set. This division was performed in such a 
manner that data coordinates of regression graph represent both 
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training and test set compounds and distributed within the whole 
descriptor space of the entire dataset. Each data point of the test 
set showed closer match with the training set compounds. The 
regression model equation was derived by using statistical multiple 
linear regression approach (Table S3 in File SI). 

Model quality assessment and validation 

The successful QSAR model must be robust enough to make 
accurate and reliable predictions of the non-investigated or query 
set compounds, therefore the obtained QSAR model from the 
training set should be subsequently validated. The conventional 
validation strategy for QSAR model analysis, based on multiple 
linear regression, include the calculation of cross vahdated squared 
correlation coefficient (r^) for internal validation and the predictive 
squared correlation coefEcient (r^p„d) for external validation. Here 
in this case the r^ was 0.878, and r^p^ed was 0.5865, which 
ultimately prove the true predictability of model and the model 
was not obtained by chance only. The parameters for model 
construction were: (i) User set: AlogP, molecular weight, number 
of H-donors, number of H-acceptors, number of rotatable bonds, 
number of rings, number of aromatic rings, molecular fractional 
polar surface area, (ii) Initial model form: Least-Squares, (rii) 
Number of nearest neighbours: 20 and (iv) Dynamic smoothing 
factor: 0.5 Graphical plot between experimental and predicted 
activities (ICso |J,M) of the training and test set compounds are 
represented in Figwe 1. 

QSAR analysis: Multiple linear regression equation. By 
considering the top eight descriptors following regression equation 
derived: 



PredictedIC5o(nM) = - 14.38 

— 3.084 xALogP 

+ 0.06609 X Molecular weight 

— 2.684 X Number of H — donors 

— 7.081 X Number of H — acceptors 
+ 8.984 X Number of rotatable bonds 

— 27 . 5 5 X Number of rings 

+ 39.81 X Number of aromatic rings 

+ 39.81 X Molecular fractional polar surface area 



Regression coefficient(N= 124,r^ =0.878,r^ 



(predicted) 



= 0.5865) 



It is evident from the above mentioned equation that the 
molecular descriptors, ALogP, number of H-donors, number of H- 
acceptors and number of rings showed negative correlation with 
respect to the biological activity. On the other hand molecular 
weight, number of rotatable bonds, number of aromatic rings and 
molecular fractional polar surface area showed positive correlation 
with the biological activity. 

Virtually designed derivatives of Capsazepine 

Forty-four capsazepine derivatives were virtually designed and 
their probable activities were predicted by the developed QSAR 
model. Virtual screening through derived QSAR model resulted 
six best hits e.g., CPZ-29, CPZ-30, CPZ-33, CPZ-34, CPZ-35 and 
CPZ-36. Predicted IC^q (l-tM) of capsazpine derivatives are 
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CPZ-36 



Figure 2. The chemical structure of capsazepine and its active derivatives namely, CPZ-29, CPZ-30, CPZ-33, CPZ-34, CPZ-35 and 
CPZ-36. 

doi:l 0.1 371/journal.pone.01 00797.g002 



summarized in Table 1 and structure of active capsazepine 
derivatives are showed in Figure 2. 

Ethics Statement 

Primary macropliage cells were isolated from the peritoneal 
cavities of the healthy female Swiss albino mice as per the 
approved protocol (AH-20 13-06) by the Institutional Animal 
Ethics Committee (lAEC) of Central Institute of Medicinal and 
Aromatic Plants, Lucknow followed by the Committee for the 
Purpose of Control and Supervision of Experimental Animals 
(CPCSEA), New D("lhi, Government of India (Registration No: 
400/0 1/AB/CPCSEA). 

Primary cell culture and treatment 

Primary cell culture was carried out as described previously 
[25]. In brief, the macrophage cells were collected from the 



peritoneal cavities of mice (8-week-old female Swiss albino mice) 
after an intra-peritoneal (i.p.) injection of 1.0 mL of 1% peptone 
(BD Biosciences, USA) 3 days before harvesting. Mice were 
euthanized by cervical dislocation under ether anesthesia and 
peritoneal macrophages were obtained by intra-peritoneal injec- 
tion of Phosphate Buffer Saline (PBS), pH-7.4. Membrane debris 
was removed by filtering the cell suspensions through sterile gauze. 
The viability of the cells was determined by trypan blue exclusion 
and the viable macrophage cells at the concentration of 0.5x106 
live cells/mL were used for the experimentation. The cells were 
suspended in RPMI 1640 medium (Sigma-Aldrich, USA) 
containing 10% heat-inactivated fetal calf serum (Gibco, USA), 
100 Hg/mL of penicillin and 100 (tg/mL of streptomycin and 
incubated in a culture plate (Nunc, Germany) at 37°C in 5% CO2 
in an incubator. Non-adherent cells were removed after 4 h by 
removing the culture media and the adherent cells were re- 
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Figure 3. Effect of capsazepine on production of TNF-a in LPS-induced inflammation in macropKiage cells (n = 5; p<0.05; * Vehicle 
versus Treatment). 

doi:1 0.1 371 /journal.pone.01 00797.g003 



suspended in RPMI 1640 medium containing 10% heat- inacti- 
vated fetal calf serum. Cells were pretreated with 1, 2.5, 5 and 
10 |J,g/mL of test compounds and standard anti- inflammatory 
drug, Dexamethasone (Sigma Aldrich, USA) at 5 |J.g/mL for 
30 min. The cells were stimulated with lipopolysaccharide (LPS, 
0.5 |J-g/mL). After incubation with LPS for 24 h, supernatants 
were collected and immediately frozen at — 80°C. Harvested 
supernatants were tested for quantification of pro-inflammatory 
mediator TNF-a by ELISA method according to the manufac- 
turer's instructions (BD Biosciences, USA). 

Quantification of pro-inflammatory cytokines 

Quantification of TNF-ot at protein level in cell culture 
supernatant was carried out using Enzyme Immuno Assay (EL\) 
kits from BD Biosciences, USA following the manufacturer's 
protocol. Briefly, the ELISA plates (96 well) were coated (100 |J,L 
per well) with specific mouse TNF-ot capture antibody respectively 
and incubate overnight at 4°C. The plate was blocked with 
200 |xL/well assay diluents. Cell Culture supernatant and standard 
(100 [iL) were added into the appropriate coated wells and 



incubated for 2 h at room temperature (20 — 25°C). After 
incubation, the plates were washed thoroughly 5 times with wash 
buffer. 100 |iL of detecting solution (detection antibody and 
streptavidin HRP) was added in to each well. Seal plate and 
incubate for 1 h at RT and then the plates were washed 
thoroughly 5 times with wash buffer. 100 |J,L of tetramethyl 
benzidine (TMB) substrate solution to each well and incubate plate 
(without plate sealer) for 30 min at room temperature in the dark. 
Add 50 |J,L of stop solution (2N H2SO4) to each well. The color 
density was measured at 450 and 570 nm using a microplate 
reader (Molecular Devices, USA). Subtract absorbance at 570 nm 
from absorbance 450 nm. The values of TNF-ot were expressed as 
|Xg/mL and the IC50 values was calculated from vector defined by 
percentage inhibition values obtained against concentration 
gradient ranging from 1-10 |J,g/mL Representative results are 
depicted in Figure 3 and Table 2. 

In vitro cell cytotoxicity 

Cytotoxicity was carried out in macrophage cells using 3-(4,5- 
dimethylthiazol-2-yl)-2,5-diphenyltetrazolium (MTT) assay. Peri- 



Table 2. Experimental in vitro activity for capsazepine against TNF-ot. 





Compound 


LPS (0.5 |ig/mL) 


Concentration (fig/mL) 


TNF-a (pg/mL) Mean+SEM 


Percentage Inhibition 


IC50 (tilVl) 


Normal 






24.60 ±4.62 






Vehicle 






1107.33±14.18 


0% 




Capsazepine 




1 HQ 


1028.47± 18.81 


7.1% 


26.7982** 




J 


2.5 (ig 


774.67 ±58.06* 


30% 








5 ng 


642.47 ±18.56* 


41 .9% 






J 


10 tig 


566.40 ±30.83* 


48.9% 




Dexamethasone 


1 


5 fig 


167.87±8.42* 


84.84% 





Note: *indicate significant effect, 

**IC5o values was calculated from vector defined by percentage inhibition values obtained against concentration gradient ranging from 1-10 ^g/mL. 
doi:1 0.1 371/journal.pone.01 00797.t002 
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Figure 4. Effect of capsazepine on percent (%) cell viability of macrophage cells using MTT assay. 

doi:10.1371/journal.pone.0100797.g004 



toneal macrophage cells (0.5x10*' cell/well) isolated from mice 
were suspended in RPMI 1640 medium (Sigma, USA) containing 
10% heat-inactivated fetal bovine serum (Gibco, USA) and 
incubated in a culture 96 well plate at 37°C in 5% COj in an 
incubator and left overnight to attach. Cells treated with 1 % 
DMSO served as a vehicle control for cell cytotoxicity study. 
Extracts and capsazepine were dissolved in DMSO. Cells were 
treated (1, 2.5, 5, 10 Hg/mL) and incubated for 24 h at 37°C in 
5% CO2. After incubation cells with treatment, 20 jxl aliquots of 
MTT solution (5 mg/ mL in PBS) were added to each well and left 
for 4 h. Then, the MTT containing medium carefully removed 
and the cells were solubUised in DMSO (100 |a,L) for 10 min. The 
culture plate was placed on a micro-plate reader (Spectromax plus 
384 with Softmax pro v5.3 software; Molecular device, USA) and 
the absorbance was measured at 550 nm. The amount of color 
produced is directly proportional to the number of viable cells. 
Cell cytotoxicity was calculated as the percentage of MTT 
absorption as follows: Percentage (%) of survival = (mean 
experimental absorbance/mean control absorbance xlOO) [26]. 
The MTT assay results are summarized in Figure 4 and Table 3. 



Table 3. MTF assay for capsazepine toxicity evaluation. 



Statistical analysis 

Results were presented as the means ± SEM and analyzed using 
GraphPad Prism 4. The ANOVA followed by turkeys multiple 
comparison tests was used to assess the statistical significance of 
vehicle verses treatment groups. Results are presented as the 
means ±SEM. Differences with a p value <0.05 were considered 
significant. IC50 values were calculated from vector defined by 
percentage inhibition values obtained against a concentration 
gradient ranging from 1 — 10 |J,g/ mL. 

IVlolecular docking 

The docking study of selected target and ligands was done by 
using Autodock Vina vO.8 (Molecular Graphics Lab at The 
Scripps Research Institute, La JoUa, CA 92037, USA). The 3D 
crystallographic structure of anti-inflammatory protein target 
tumor necrosis factor-a (TNF-ot) was retrieved through Brook- 
haven Protein DataBank (PDB) (http://www.pdb.org) (PDB ID: 
2AZ5). The crystallographic protein structure of TNF-a complex- 
es with known inhibitor was selected for docking procedure 
validation by re-docking approach and also to know the standard 
docking energy and binding site. The valency and hydrogen bonds 
of the ligands as well as target protein was subsequently satisfied. 



Compound 


LPS (0.5 ^g/mL) 


Concentration (^g/mL) 


% Live Cells Mean+SEM 


Normal 






100.00±4.53 


Vehicle 


J 




105.89±2.20 


Capsazepine 


J 


1 


103.53±2.92 




J 


2.5 


99.07±2.32 




J 


5 


93.87 ±2.46 




J 


10 


89.47 ±1.92 


Dexamethasone 


J 


5 


102.08±2.02 



Note: Significant (p<0.05) changes in % live cells were not observed in any treatment. 
doi:l 0.1 371/journal.pone.0100797.t003 
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Table 4. Comparison of binding affinity of capsazepine and its active derivatives in terms of docl<ing energy and binding site 
residues against anti-inflammatory receptor TNF-a (PDB: 2AZ5). 





Compounds 


Docking binding 
energy (l<cal/mol) 


Chain ID 
(dimer form) 


Binding pocliet residues 
within 4 A radius 


interacting residues and 
length (4 A) 


Re-docking of bound inhibitor-307* 
in TNF-Q( dimer complex form 


-9.1 


Chain A: 


LEU-57, TYR-59, SER-60, TYR-119, 
LEU-120, GLY-121, GLY-122, TYR-151 








Chain B: 


LEU-57, TYR-59, TYR-1 19, LEU-120, 
GLY-121 




Capsazepine 


-7.3 


Chain A: 


LEU-59, TYR-119, LEU-120, GLY-121 








Chain B: 


LEU-57, TYR-59, SER-60, GLN-61, 
TYR-119, LEU-120, GLY-121, TYR-151 


Chain B: LEU-120 (2.3) 


CPZ-29 


-7.2 


Chain A: 


LEU-57, TYR-59, TYR-1 19, LEU-120, 
GLY-121, GLY-122, TYR-151, ILE-155 








Chain B: 


LEU-57, TYR-59, SER-60, TYR-119, 
LEU-120, GLY-121, TYR-151 




CPZ-30 


-7.9 


Chain A: 


TYR-59, SER-60, GLN-61, TYR-119, 
LEU-120, GLY-121, TYR-151 








Chain B: 


TYR-59, SER-60, TYR-119, LEU-120, 
GLY-121, TYR-151 


Chain B: LEU-120 (2.3) 


CPZ-33 


-6.8 


Chain A: 


TYR-59, SER-60, TYR-119, LEU-120, 
GLY-121, TYR-151 


J 






Chain B: 


TYR-59, SER-60, TYR-119, LEU-120, 
GLY-121 


Chain B: LEU-120 (2.7) 


CPZ-34 


-7.8 


Chain A: 


TYR-59, TYR-119, LEU-120, GLY-121 


Chain A: LEU-120 (2.7) and 






Chain B: 


LEU-57, TYR-59, SER-60, GLN-61, 
TYR-119, LEU-120, GLY-121, TYR-151 


Chain B: TYR-151 (3.2) 


CPZ-35 


-7.3 


Chain A: 


LEU-57, TYR-59, SER-60, TYR-119, 
LEU-120, GLY-121TYR-151 








Chain B: 


LEU-57, TYR-59, SER-60, TYR-119, 
LEU-120, GLY-121, TYR-151 




CPZ-36 


-6.9 


Chain A: 


TYR-59, SER-60, TYR-119, LEU-120, 
GLY-121, TYR-151 








Chain B: 


TYR-59, SER-60, TYR-119, LEU-120, 
GLY-121, TYR-151 





Note: "-" represent no H-bond and • refer TNF-a dimer inhibitor-307 name (6,7-DIMETHYL-3-[(METHYL{2-[METHYL({l-[3-(TRIFLUOROMETHYL)PHENYL]- lH-INDOL-3- 
YL}METHYL)AMINO]ETHYL}AMINO)METHYL]- 4H-CHROMEN-4-ONE). 
doi:l 0.1 371 /journal.pone.Ol 00797.t004 



An extended PDB format, termed as PDBQT file was used for 
coordinate files that includes atomic partial charges. The software 
automatically convert PDB file into PDBQT that was further used 
for docking [27]. Polar hydrogen atoms were added to the protein 
target to achieve the correct ionisation and tautomeric states of 
amino acid residues such as HIS, ASP, SER and GLU. For 
software standardization, native ligand of co-crystallised complex 
was first extracted and re-docked to its corresponding binding site 
using AutoDock Vina vO.8. The docking of ligand with receptor 
TNF-ot was performed for three times and the average of 
consecutive results was taken as final binding/ docking energy. 
Binding pose with the lowest docked energy belongs to the top- 
ranked cluster was selected as the final model for post-docking 
analysis with UCSF Chimera vl.5.3 (NCRR, University of 
California, San Francisco, supported by NIH P41 RR001081) 
and PyMOL (The PyMOL Molecular Graphics System, Version 
1.5.0.4 Schrodinger, LLC, USA). 

Screening through pharmacokinetic properties 

Most of drugs in development failed during clinical trials due to 
poor pharmacokinetics parameters [28,29] . These properties such 
as absorption, distribution, metabolism, excretion and toxicity 
(ADMET) are important in order to determine the success of the 



compound for human therapeutic use. Some important chemical 
descriptors correlate well with ADMET properties such as polar 
surface area (PSA) as a primary determinant of fraction 
absorption, low molecular weight (MW) for oral absorption. The 
distribution of the compound in the human body depends on 
factors such as blood-brain barrier (Log BB), permeability such as 
apparent Caco-2 permeability, apparent MDCK cell permeability. 
Log Kp for skin permeability, volume of distribution and plasma 
protein binding (Log Khsa for Serum protein binding). It has been 
reported that excretion process that eliminates the compound from 
human body depends on the molecular weight and octanol-water 
partition coefficient (LogP). Similarly, rapid renal clearance is 
associated with small and hydrophUic compounds. The metabo- 
lism of most drugs that takes place in the liver is associated with 
large and hydrophobic compounds. Higher lipophUicity of 
compounds leads to increased metabolism and poor absorption, 
along with an increased probability of binding to undesired 
hydrophobic macromolecules, thereby increasing the potential for 
toxicity. In spite of some observed exceptions to Lipinski's rule, the 
property values of the vast majority (90%) of the orally active 
compounds are within their cut-off limits [30 — 31]. In addition, 
the bioavailability of derivatives was assessed through topological 
polar surface area analysis. We calculated the polar surface area 
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Figure 5. Docking pose of capsazepine and its active derivatives on anti-inflammatory receptor TNF-a (PDB: 2AZ5). (a) Docking 
protocol standardization by re-docking of co-crystallized ligand on TNF-a with docking energy —9.2 kcal/mol, (b) Capsazepine docked on TNF-oc with 
binding energy —7.3 kcal/mol, (c) CPZ-29 docked on TNF-a with binding energy —7.2 kcal/mol, (d) CPZ-30 docked on TNF-a with binding energy - 
7.9 kcal/mol, (e) CPZ-33 docked on TNF-a with binding energy -6.8 kcal/mol, (f) CPZ-34 docked on TNF-ot with binding energy —7.8 kcal/mol, (g) 
CPZ-35 docked on TNF-a with binding energy —7.3 kcal/mol and (h) CPZ-36 docked on TNF-a with binding energy —6.9 kcal/mol. 
doi:1 0.1 371 /journal.pone.01 00797.g005 



(PSA) by using method based on the summation of tabulated 
surface contributions of polar fragments termed as topological 
PSA (TPS A). Generally, it has been seen that passively absorbed 
compounds with a PSA> 1 40 A" are thought to have low oral 
bioavailability. Calculations of other important ADME properties 
of capsazepine derivatives were performed through Discovery 
Studio v3.5, USA (2013). We also screened capsazepine and its 
derivatives through TOPKAT toxicity estimation using Discovery 
Studio v3.5. TOPKAT computes a probable value of toxicity for a 
submitted chemical structure from a quantitative structure-toxicity 



relationship (QSTR) equation. The product of a structure 
descriptors and its corresponding coefficient is the descriptors 
contribution to the probable toxicity. 

Results and Discussion 

In the present work, derivatives of capsazepine were evaluated 
for their anti-inflammatory activity through the developed QSAR 
model and docking studies. Structure activity relationship has been 
denoted by QSAR model showing significant internal prediction 
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Table 5. Compliance of capsazepine and its derivatives to the computational parameters of oral bioavailability (Lipinski's rule of 
five). 



Compound 


Oral l^i/^nw stilnl^ili'l-if' TDCA 
widl UIUdvclllcli.ll Illy • I r9n 

«140 A^) 


MW «S00) 


ALog P (£5) 


n uuiiu uuilur 
(S5) 


ri uuiiu ci(.i.c|jiur 
(S10) 


Rule of 5 violations 


Capsazepine 


87.82 


376.909 


5.393 


3 


3 


0 


CPZ-29 


74.47 


489.113 


8.501 


1 


3 


1 


CPZ-30 


74.47 


433.007 


6.542 


1 


3 


1 


CPZ-33 


83.7 


400.534 


5.163 


1 


4 


0 


CPZ-34 


116.7 


358.455 


4.486 


4 


4 


0 


CPZ-35 


74.47 


517.166 


9.414 


1 


3 


2 


CPZ-36 


74.47 


461.06 


7.589 


1 


3 


1 



Footnote: A, absorption; D, distribution; M, metabolism; E, excretion; IPSA, topological polar surface area; MW, molecular weight; LogP 
coefficient, a measure for lipophilicity. 
doi:l 0.1 371/journal.pone.Ol 00797.t005 



octanol/water partition 



(r^ = 0.8779) and external prediction (r^p^ed = 0.5865) (Figure 1). A 
total of 124 known inhibitors of TNF-ot were used for QSAR 
modeling against 57 2D and 3D chemical descriptors. Eight 
descriptors were found to be significantly responsible for anti- 
inflammatory activity (Table 1). A forward feed multiple linear 
regression QSAR model was developed. Both internal and 
external validations were performed for the developed model. A 
low residual value for each compound in the dataset defines the 
degree of correlation between observed and predicted values and 
the models predictive ability. Screening of derivatives through 
developed QSAR model indicated that derivatives CPZ-29, CPZ- 
30, CPZ-33, CPZ-34, CPZ-35 and CPZ-36 showed significant 
activity in compared to capsazepine's in vitro IC50 against TNF-a 
(Table 2). In 1989, it was found that capsazepine and some of its 
derivatives possessed extraordinary anti-inflammatory activity 
against COX-2 and inducible nitric oxide (iNO) [11,12]. In the 
present work, we report anti-inflammatory activity of 44 virtually 
designed capsazepine derivatives with lactone ring pharmaco- 
phore against pro-inflammatory target TNF-ot. The activity of 
newly designed derivatives were predicted through the developed 
QSAR model and the derivatives CPZ-33 and CPZ-34 found to 
be better in activity as compared to capsazepine, whereas CPZ-30 
was found close to capsazepine. AH six derivatives and parent 
compound (Figure 2) were further selected for in silico target- 
receptor interaction and ADMET studies. 



Binding affinity study through docking against TNF-ot 

The aim of the molecular docking study was to elucidate 
whether capsazepine derivatives modulate the anti-inflammatory 
target and also to identify the binding site against well-known 
human anti-inflammatory molecular target TNF-a. The native 
ligand re-docking study indicates that the software predict the 
reliable results. Thereafter the predicted active derivatives were 
subjected to molecular docking studies. The docking results 
provided pertinent information about the binding affinity, binding 
energy and orientation of ligand-receptor interactions. The 
docking results are summarized in Table 4. It has been found 
that capsazepine and its active derivatives bound to the same 
active site as reported in the PDB protein crystallographic 
structure database. Capsazepine showed significant binding 
affmity to TNF-ot dimeric structural unit (A and B chain residues) 
with binding energy of —7.3 kcal/mol, similarly, capsazepine 
active derivatives CPZ-34, and CPZ-30 showed high binding 
affmity to TNF-ot dimeric structural unit with docking energy of — 
7.8 and —7.9 kcal/mol, respectively. Docking pose of capsazepine 
and its active derivatives on receptor TNF-ot are showed in 
Figure 5. On the other hand, capsazepine derivatives CPZ-29 and 
CPZ-35 showed moderate binding affinity to TNF-ot dimeric 
structural unit with binding energy of —7.2 and —7.3 kcal/mol, 
respectively. The capsazepine derivatives CPZ-33 and CPZ-36 
showed low binding affmity to TNF-ot dimeric structural unit with 
binding energy of —6.8 and —6.9 kcal/mol, respectively. The 



Table 6. Compliance of capsazepine and its active derivatives to the computational parameters of pharmacokinetics {ADME). 



Aqueous BBB Intestinal 
Compound AlogP PSA_2D Plasma protein binding Hepatotoxicity CYP2D6 binding solubility penetration absorption 



Capsazepine 


4.674 


57.793 


True (highly bounded) 


True (toxic) 


False (non-inhibitor) 


2 (low) 


1 (good) 


0 (good) 


CPZ-29 


7.929 


34.022 


True (highly bounded) 


False (non-toxic) 


False (non-inhibitor) 


1 (poor) 


4 (undefined) 


3 (very poor) 


CPZ-30 


5.969 


34.022 


True (highly bounded) 


False (non-toxic) 


False (non-inhibitor) 


1 (poor) 


0 (very good) 


1 (moderate) 


CPZ-33 


4.591 


42.952 


True (highly bounded) 


True (toxic) 


False (non-inhibitor) 


2 (low) 


1 (good) 


0 (good) 


CPZ-34 


3.914 


78.609 


False (poorly bounded) 


False (non-toxic) 


False (non-inhibitor) 


3 (good) 


2 (medium) 


0 (good) 


CPZ-35 


8.841 


34.022 


True (highly bounded) 


False (non-toxic) 


False (non-inhibitor) 


0 (very poor) 


4 (undefined) 


3 (very poor 


CPZ-36 


7.016 


34.022 


True (highly bounded) 


False (non-toxic) 


False (non-inhibitor) 


1 (poor) 


4 (undefined) 


3 (very poor 



Abbreviations; AlogP, the logarithm of the partition coefficient between n-octanol and water; PSA, polar surface area, CYP450 cytochrome P450, PPB plasma protein 
binding, BBB blood brain barrier. 
doi:1 0.1 371 /journal.pone.Ol 00797.t006 
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Figure 6. Plot of polar surface area (PSA) versus ALogP for capsazepine and its derivatives showing the 95% and 99% confidence 
limit ellipses corresponding to the blood brain barrier (BBB) and intestinal absorption. 

doi:10.1371/journal.pone.0100797.g006 



chemical nature of capsazepine & its derivatives binding site 
amino acid residues on TNF-a dimeric structural unit (Chain A & 
B) were aliphatic {e.g., LEU-57, LEU- 120, GLY-121, GLY-122), 
hydroxyl group containing [e.g., SER-60), and aromatic [e.g., TYR- 
59 and TYR-119). The capsazepine, CPZ-30, CPZ-33 and CPZ- 
34 showed H-bonds with the TNF-ot dimeric unit residues, this 
suggest high structural stability and may lead to high inhibitory 
activity of capsazepine and its derivatives on TNF-ot active site. 

Bioavailability and ADME parameters screening for drug 
likeness 

The compound's good absorption or permeation through blood 
brain barrier is measure by its Log? that must be less than 5 [22 — 
23]. Results of pharmacokinetic screening revealed that capsaze- 
pine and its most active derivatives CPZ-33 and CPZ-34 followed 
the Lipinki's rule of five for oral bioavailability. However CPZ-29, 
CPZ-30, and CPZ-36 showed lipophilic nature due to high LogP 
value, while compound CPZ-35 showed both high lipophilicity 
and low membrane permeability due to high LogP and molecular 
weight. These ADMET screening results are summarized in 
Table 5 and 6. The ADME descriptors of capsazepine and its 
derivatives were calculated for drug likeness studies. The intestinal 
absorption and blood brain barrier penetration were predicted by 
developing an ADME model using descriptors 2D PSA and 
AlogP98 that include 95% and 99% confidence ellipses. These 
ellipses define regions where well-absorbed compounds are 



expected to be found. The results of DS-ADME model screening 
showed that capsazepine derivatives CPZ-33 and CPZ-34 possess 
99% confidence levels for human intestinal absorption and blood 
brain barrier (BBB) penetration. Similarly, another predicted 
active capsazepine derivative CPZ-30 also showed 99% confi- 
dence level for intestinal absorption and 95% confidence level for 
BBB penetration. The capsazepine derivatives CPZ-29, CPZ-35 
and CPZ-36 fall outside the ADME model ellipses filter, which 
indicate its poor intestinal absorption and BBB penetration ability. 
The plot of polar surface area and ALogP fpr capsazepine and its 
derivatives are represented in Figure 6. 

Toxicity risks assessment 

The USFDA (US FDA, United States Food and Drug 
Administration) standard toxicity risk predictor software TOP- 
KAT (Discovery Studio, Accelrys, USA) locates fragments within 
the compound that indicate a potential threat to toxicity risk [26] . 
Toxicity screening results of TOPKAT for capsazepine and its 
derivatives showed that studied compounds possess no risk of 
carcinogenicity, mutagenicity and skin irritation, however it 
possess high developmental or reproductive toxicity potential at 
high doses or long term therapeutic use in human. The 
capsazepine and its derivatives CPZ-30, CPZ-33 and CPZ-34 
showed strong skin sensitization capacity. Similarly, capsazepine 
derivatives CPZ-33 and CPZ-34 also showed mild ocular irritancy. 
Other detail predicted toxicity parameters are summarized in 
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Table 6. The results of toxicity risk for capsazepine and its active 
derivatives showed moderate to good drug score, in compared 
with capsazepine (Table 7a). Similarly, toxicity screening results of 
USFDA rodent carcinogenicity, Ames mutagenicity, developmen- 
tal toxicity potential, aerobic biodegradability, ocular irritancy and 
skin irritancy also showed positive response to capsazepine and its 
derivatives (Table 7b). 

Anti-inflammatory potential of Capsazepine 

To examine the effects of capsazepine on LPS-induced pro- 
inflammatory cytokine TNF-ot in macrophages, culture superna- 
tant from various treatment groups were used to determine their 
production. Treatment of capsazepine inhibited (p<0.05) the 
production of LPS-induced inflammatory mediator in dose 
dependant manner. (Figure 3, Table 2). 

Measurement of the cell viability 

The in vitro effect of capsazepine on cell viability in peritoneal 
macrophage cells isolated from mice was evaluated using MTT 
assay. The significant change in percent live cell population was 

not observed (p<0.05) at any concentration of the treatment when 
compared with normal cells (Figure 4, Table 3). 

Conclusions 

The experimental in vito evaluation of capsazepine against pro- 
inflammatory mediator TNF-a indicated that capsazepine mediate 
significant inhibitory effect on the TNF-a. The predicted activity 
of capsazepine was comparable with the experimental results. 
Ligand-based virtual screening through developed QSAR model 
resulted in six best hits for capsazepine derivatives CPZ-29, CPZ- 
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